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Abstract 

Cylindrical implosions driven by intense heavy ions beams should be instrumental in a near 
future to study High Energy Density Matter. By rotating the beam by means of a high frequency 
wobbler, it should be possible to deposit energy in the outer layers of a cylinder, compressing the 
material deposited in its core. The beam temporal profile should however generate an inevitable 
irradiation asymmetry likely to feed the Rayleigh- Taylor instability (RTI) during the implosion 
phase. In this paper, we compute the Fourier components of the target irradiation in order to 
make the junction with previous works on RTI performed in this setting. Implementing a ID 
and 2D beam models, we find these components can be expressed exactly in terms of the Fourier 
transform of the temporal beam profile. If T is the beam duration and its rotation frequency, 
"magic products" can be identified which cancel the first harmonic of the deposited density, 
resulting in an improved irradiation symmetry. 

PACS numbers: 41.75.2i, 62.50.1p, 52.59.2f 
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I. INTRODUCTION 



The study of matter under extreme conditions of pressure and density is of great inter- 
est for astrophysics, planetary sciences and inertial fusion {l-4|. Particularly appealing is 
the perspective to produce cylindrical implosions with a high degree of symmetry. Such 
experiments would allow to generate large volumes of high-energy density matter (HEDM), 
including strongly coupled plasmas. Among others, the experimental realization of the 
once predicted metallic state of Hydrogen [5] is envisioned at the Gesellschaft fiir Schw- 
erionenforschung (GSI) near Darmstadt, Germany. Within the framework of the Facility 
for Antiproton and Ion Research (FAIR) currently under construction at GSI [6], the so- 
called LAPLAS experiment (LAboratory of PLAnetary Sciences) aims at implementing such 
cylindrical scheme to study equation of state and transport properties of HEDM 

Figure 1 sketches the typical experimental scheme implemented for cylindrical implosions. 
A hollow ion beam hits the absorber part of the cylinder. By tailoring its energy so that 
the Bragg peak of the ions falls well outside the cylinder, it is possible to obtain a very 
homogenous energy deposition inside the cylinder. The following compression has been the 
object of various numerical and analytic studies in recent years, to assess, among other. 
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The hollow beam generation should be achieved by a high frequency wobbler rotating 
the ion beam [15]. The rotation frequency necessary to achieve an acceptable degree of 
irradiation symmetry has been calculated by Piriz et al [12] and found in the GHz range. 
For such a fast deposition, the target motion during the deposition can be neglected. But 
even so, an inevitable source of asymmetry comes from the beam temporal profile itself. 
Because the beam duration is finite, the number of ions deposited along the absorber will be 
inhomogeneous. The overall effect of such asymmetry on the pressure driving the compres- 
sion has been investigated in Refs. jll| and jisl, where the amplitude of the asymmetry has 
been evaluated. But previous works on the RTI for this setting have showed that stability is 
a matter of both the excited wavelength, and the excitation amplitude. It is thus necessary 
to know the Fourier components of the ion deposition in the pusher, precisely because RTI 
analysis is performed in Fourier space. This procedure is reminiscent of studies conducted 
for Inertial Confinement Fusion, where the driver energy deposition on the spherical pellet 




FIG. 1. (Color online) Typical cylindrical implosion scheme. Not to scale. 



has to be Fourier analyzed to assess its 
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The goal of this paper is to compute the Fourier components of the ion deposition in the 
absorber. We start implementing a ID beam model and express the harmonics amplitude of 
the ion density deposition, in terms of the Fourier transform of the temporal beam profile. 
We then consider a beam spatially extended in the transverse direction, and generalize the 
ID results. We then turn to the stability analysis, making the junction with previous works. 
Finally, we address the possibility of canceling the first harmonic. In the case of a parabolic 
temporal profile, an infinite series of product are found for which the first harmonic 
vanishes exactly. 

II. ID BEAM MODEL 

A. Irradiation spectrum calculation 

Consider a ID ion beam with temporal profile /(t), normalized to unity so that J I{t)dt = 
1. The beam is rotated by a wobbler over the absorber region at velocity Vt rad/s (see Figure 
[2j). The number of ions brought by the beam between t and t + rft is I{t)dt. Suppose ions 
arriving at t = hit the target at ^ = 0. Then those arriving t = 6 /Vt hit the target at 
angle 0, as explained on the Figure. The number of ions deposited between 9 and 9 + d9^ 
from time t to time t + dt is therefore. 



Then, ions arriving 21t{ /Vt seconds later, or before, hit the very same point, for any / integer, 
positive or negative. The total number of ions eventually deposited between 9 and 9 -\- d9 is 
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FIG. 2. (Color online) The beam intensity profile is I{t). It rotates at velocity rad/s along the 
circle. Ions arriving at t = + 2k7v/ft all hit the same point of the circle. 
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As expected, the density deposition dN{6)/d6 = p{6) is periodic, of period 27r. 

We now turn to the Fourier transform of the density deposition p(6) in the absorber ring. 
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Having a periodic function, it could be possible to Fourier analyze it integrating from to 
27r. However, integrating from — oo to oo is more general and allows to deal with the Fourier 
Transform of the beam profile as well. We may now permute the infinite sum with the 
integral because we are dealing with "real world" regular functions for which such operation 
can be done. Making the substitution u = 9/Q + 21ti /Vt^ we find 
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The term between parenthesis is just the Fourier transform of the beam temporal profile for 
the frequency value "sf^" . We thus have 

CO oo 
/ = — OO l = — oo 

According to some properties of the "Dirac's comb" function ^^(^(s — /) [3], we have 

oo oo 

/=— oo /=— oo 

and we finally come to, 

oo oo 

p{s) = E T{sn)e-'^^-^ = T{sn) E Sis - I). (7) 

/=— OO /=— oo 

The infinite sum "^^(8 — 1) functions at the right hand side implies the spectrum is discrete. 
This comes from the fact that we consider the Fourier transform of the periodic function 
p{9). Furthermore, the sum of 5's implies that the discrete values of the spectrum are the 
integers / = — oo . . . oo. The amplitudes of the harmonics are therefore, 

m = m, 

= T{n) + T{-n), 

p(2) = T{2Q) + T{-2Q), 

p{i) = T{m) + T{-m), yi e n. (s) 

From the harmonics above, the density deposition p(9) is expressed as 
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Such is the final result, where the amplitude of the harmonics of the density deposition is 
expressed in terms of the beam temporal profile. Let us now illustrate the result considering 
the example of a parabolic temporal profile. 




FIG. 3. Temporal profile of the beam from Eq. (|TT1) with T = 2 s. The number of ions hitting the 
target from t to t + dt is I{t)dt. J I{t)dt = 1. 

B. Parabolic beam temporal profile 

For a parabolic temporal beam profile normalized to 1 (see Fig. [3j), 

2" 



I{t) = — 
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for -T/2<t< T/2 



= otherwise, (11) 

we find for the n*'^ harmonic of the density deposition p{0), 

N ?■/ -12cos(n™2) 24sin(nTO/2) 

Figure H] shows the density p{6) profile calculated from Eq. ([2j) and the terms / = of Eq. 
([9j) together with sum ([9]) up to / = ±1 and ±4. The term / = is just the mean value of 
the density deposition. One can notice how the sum up to / = ±4 is already very close to 
the exact density. 

Let R be the radius of the circle where the ion beam are deposited, the n*^ harmonic 
pertains to the wavelength A = 2TiR/n and the wave vector 

Let us assume the wobbler rotates times in T seconds, that is TVt = 2tiN ^ the relative 
amplitude of the first harmonic is 

m 3 



p(o) T^^m 



0.3% for N = 10. (14) 
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FIG. 4. (Color online) Black: Density profile calculated from Eq. ([T]). Blue: Term / = in Eq. 
(j9j). Green: Eq. (ED up to / = ±1. Red: Eq. (|9j) up to / = ±4. Parameters are T = 2 s and = 5 
rad/s. 



C. General analysis for finite duration beams 



We now derive some general result regarding functions I{t) which vanish out of a finite 
time range [0,T]. According to Eq. ([8]), the n}^ harmonic amplitude of the ion density 
deposition is given 

POO pT 

p{n) = T{nn) + T{-nn) = 2 I{t) cos{nnt)dt = 2 I{t) cos{nnt)dt. (15) 



Integrating by part, we find 

2 2 

= -7^ \Ht) sin(nm)]^ / r(t) 8m(nnt)dt, (16) 

nil nil Jq 

where r{t) = dl/dt. Repeating again the process gives, 

Pin) = ^ [lit) sin(nm)]J^ + [I' it) cosinQt)]^ ~ ^"^^^ cos{n^t)dt. (17) 

It is thus straightforward that if /(O) = I{T) = 0, while the first derivative does not vanish 
in t = or T, p{n) behaves like 1/n^. More generally, p{n) behaves like if all the 

derivative up to the k^^ vanish in t = and T. We thus come to the same conclusion than 

n 

Ref. pjj|. Additionally, we find that harmonic amplitudes can be tuned by choosing QT 
appropriately. This important point is detailed in Section |Vl 



III. 2D BEAM MODEL 



A. Irradiation spectrum calculation 

We have so far considered a point-like beam with a temporal profile. We thus need to 
account for a beam with spatial extension. The beam transverse density is defined in its 
own rest frame by 

/(t,u)=^(t)a(u) = e(t)(T(«), (18) 

where u is the distance from the beam center to a given beam point. We thus consider 
beams with cylindrical symmetry. 

Let us calculate the total number of ions arriving at the target point localized by (i?, 9) 
(the red point on Fig. [5]). We can consider here that the target is hit by various beamlets. 
The target point (i?, 9) is hit by n "beamlets" represented by the blue circles. Of course, 
n is eventually infinite. We thus start the analysis with n finite, before we make it tend to 
infinity. 

Reasoning first in the target polar coordinates (the green elements on Fig. Ej), the number 
of ions brought by the first beamlet is, 

dN = a{ui)uid9u du (19) 

By analogy with Eq. (Ej) from the ID calculation, we deduce the ions density deposited by 
this beamlet, 

dpi{R,9) = ^ ^ ^ ^ (y{ui)uid9udu. (20) 

All the beamlets pictured on Fig. [5] up to "beamlet n" will deposit ions in {R^9). But it 
is clear that beamlets with n > 1 deposit their ions at target point (i?, 9) some time later, 
depending on how large is the angle a the beam needs to rotate for those beamlets to cover 
the target point. Therefore, the beamlet number fc, deposits 

dPk{R. ^) = ^Y.^[ ) ^(^)*^^- (21) 

l=-oo ^ ^ 

Summing the contributions from all beamlets, we integrate over the arc of radius R inter- 
secting the beam spot. It is of course much more convenient to work in polar coordinates 
with respect to the target center instead of the beam spot center. We thus integrate the 
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FIG. 5. (Color online) Calculation for a beam with temporal and transverse spatial profile. Not 
to scale. 

expression above from a = to 9^- When switching from {u^Ou) to {R^a) coordinates, 
the differential element dOudu changes to 6)drd6^ where J is the Jacobian of the trans- 
formation. We thus find 



a=-dn 



C ( ^ ^ ) ^ [^(-^' ")] ^)JiR^ dr. (22) 



This result simplifies substantially due to the identity proven in Appendix [A] 



u{R^ a)J{R^ a) = i?, 



(23) 



SO that 



p{R,e) 



Rdr 
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a=-0r7 



^ [ ) a [u{R, a)\ da. 



l=—oo 



(24) 



Note that 9)/ Rdr can be interpreted as the ion density deposited in the absorber along 
the circle on radius R. Also, integration from a = —9^ to 9^ can be replaced by an 
integration from —oo to oo since the integrand vanishes out of [— ^rn, ^m]- Om interest now 
lies in the Fourier transform with respect to 9 of the function above. Noteworthily, Eq. ( !24l) 
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can be written as, 



with 



if*9m 



f{9 + a)g{a)da^ 



(25) 



9{r) = a[u{R,T)] . 



r + 2/7r 



(26) 



Changing a —a, and noting that g is an even function of a, we rewrite Eq. (I25l) . 

/oo 
f{6 - a)g{a)da, (27) 
-OO 

which appears to be the convolution product of the two functions. The first of these function 
is very reminiscent of Eq. ([2]) from the ID case. It Fourier transform with respect to 9 can 
be calculated following the very same line and is found adapting Eq. ([7j) as 



(28) 



Because the Fourier transform of a convolution product is the product of the Fourier trans- 
forms, we can now write directly 



Rdr 



= P2d{R, s) 



a [u{R^ s)] 



(29) 



The interpretation of the equation above is simple: the n^^ harmonic with respect to 9 is 
now the n^^ harmonic of the point-like beam, times the n^^ harmonic of a beam form factor. 
For a point-like beam, this factor is a 5 function which Fourier transform is 1. This equation 
also shows that one should expect greater homogeneity with a spread out beam since both 
harmonics (temporal and spatial) are multiplied to yield the deposition harmonics. 



B. Parabolical beam in time and space 



We now apply the formalism developed previously to a parabolical beam both in time 
and space. The temporal profile ^{t) is thus identical to the function defined by Eq. f fTTI) . 
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and its harmonics are 



i{s) = 12 



2sin(sT/2) - sTcos(sT/2) 



The spatial profile a{u) is 



Rb 



otherwise. 



for — Rh < u < Rh 



Choosing here = (Rc, R), the vectorial relation Rc + u = R allows to derive, 

i{R,9) = ^R^ - Rl-2RcUcose, 



and, 



a[u{R,e)\ = 2 



Rl- R^ - Rl + 2rRcCose 



nRt 



with r G [Rc — Rh'>Rc + Rh\ and 9 G [— ^m] with 



cos 6 J, 



it -\- It^ — It^ 
2RRr 



The Fourier transform of a[u{R^ 9)] is 

/oo 
-oo 

The result reads, 



RI a[u{x, s)] = -4- 



^^^^2x^(i+^^)^|M£!^ ^^g^^Q^ + (e2 - 1 - x2) sin(se) 



where, 



a; = — e[l — e, 1 + eJ, e = — , B = arccos 
itc Rc 



1 + - 

2a^ 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



Each circle of radius r thus yields a different Fourier transform. Equation (1371) is plotted on 
Figure [6] for e = 1/5. The spectrum is continuous in s and of infinite extension because the 
source function is finite only over a finite 9 interval. Because ion deposition tends to zero 
near the borders of the beam, the harmonics amplitude vanish for x = 1 ± e. 



C. ID approximation 



Upon which condition shall the Fourier transform in the 2D model be comparable to the 
ID one, for any radius? If the form factor a[u{R^ s)] remains constant over the first harmonics 
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Outer beam border 



FIG. 6. (Color online) Plot of Equation ([37|) for e = 1/5. The inner border of the beam spot is 
defined by x = 0.8 = 1 — e, and the outer border by x = 1.2 = 1 + e. 

of the temporal profile. Indeed, these first harmonics are the most important, because their 
amplitude is larger. By property of the Fourier transform, the function a[u{R^ s)] has an 
extension As such that AsA9 ^ 1, where A9 is the extension of a[u{R^9)]. It is clear 
that maximal extension is achieved for R ^ R^^ with ^ R^j/R^. Hence, the Fourier 
transform a[t^(r, 5)] has a minimal extension ^ R^/Rb- Confirmation if found on Fig. [6l 
where a[u{x^ s)] is all the more peaked than x is close to 1, i.e. R = Rc. 

Therefore, if a[u{R^s)] can be considered constant for s < Rc/Rb while the ID model 
has its harmonic for 5 = 1,2,..., then the beam can be considered point-like for harmonics 
5 = 1... Rc/ Rh^ implying only Rc > Rb- With Rc = 5Rb for example, harmonics up to 5 
should be close to the ID model even at i? = Rc. 

Turning now to the stability analysis, we use the ID model, precisely because the most 
relevant harmonics are the first ones. Indeed, we will check that the analysis can be con- 
ducted in terms of the s = 1 harmonic only. For the ID model to be valid, we thus simply 
need Rb < Rc 



IV. STABILITY ANALYSIS 

The RTI at the absorber-pusher interface has been studied in Ref. jl^. Due to the 
expected experimental conditions, both sides of the interface are occupied by solids where 
elastic-plastic effects have to be accounted for. In this respect, it has been found that 
the RTI for a perturbation of wavelength A and initial amplitude depends on the two 
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dimensionless parameters, 



A 



47rG" 



where g is the acceleration of the interface, p, G and Y the density, the shear modulus and 
the yield strength of the pusher respectively. 

The stability diagram appears on Fig. [71 In the present problem, the excited wave- 
lengths are = 27ri?/n, n G N. The A's are thus decreasing with n, the largest one being 
Ai = 27rR. For this mode to be stable, i.e. to lie in the shaded area on Fig. [71 it is necessary 
(but not sufficient) that 

In the stability diagram, the excited modes yield points (A^,(^^). On the one hand, it is 
clear that A^ < Ai forn 7^ 1. On the other hand, harmonics amplitude are usually decreasing 
functions of n. This can be checked on Eq. (fT2l) for the parabolic time profile beam, or more 
generally on Eq. ( fT71) . To be more accurate, we can say that the harmonics as a function of 
n are contained within an envelop decreasing with n. 

As a consequence, if the point (Ai,^i) lies in the stability region, the points A^ < Ai for 
n 7^ 1 typically lie inside as well, as a result of the very shape of the shaded area. 

The equation of the upper limit of the stability region reads jl^ . 

^=1-Vx. (41) 
We thus need ^1 < 1 — y Ai for the system to be stable, which reads. 



Following Ref. [145, introduce, 

Po = pgh, (43) 



allowing to rewrite Eq. f l42l ) as 



< 1 - a//^^- (44) 



^/3Y h V ^ttG h 

We can consider here ii/h = p(l)/p(0). Using Eq. (|T^ for the parabohc time profile 

beam and Ai = 27ri?, we find 
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If 7rpoR/2Gh <C 1, which seems to be the case in reahstic situations, we just have, 



which is the criterion obtained in 



V. CANCELING THE FIRST HARMONIC 

An important issue for the present problem is the presence of the first harmonic. One 
reason is that its amphtude is usuaUy the largest one. But most of all, the RTI analysis 
we used has been developed for a planar interface. If the interface is circular, the planar 
analysis is still valid providing the wavelength of the perturbation is much smaller than the 
circumference of the circle. In the present case, the first excited wavelength Ai is precisely 
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FIG. 8. Plots of functions x and tanx vs. x = The equation x — tanx has an infinite 

number of solutions. 

the circumference. It is thus probable that the planar analysis fails for this mode. Canceling 
the first harmonic could solve the two problems at once: the planar analysis would just be 
applied to the next modes Ai/2, Ai/4 . . with a better accuracy and also, we could reduce 
significantly the irradiation asymmetry. 

From Eq. f fTSl ). we see than canceling the first harmonic means tailoring the pulse shape 
I{t) and the rotation frequency Vt to fulfill, 

rT/2 



I{t) co8{Qt)dt = 0, 

T/2 

where we choose t = as the middle of the pulse. Substituting u = Qt gives 



(47) 



QT/2 



I{t/n) cos{u)du = 0. 



(48) 



Since / is necessarily positive, the condition cannot be fulfilled if the cosine function takes 
only positive values in the integration range, i.e. QT/2 < 7r/2. We thus have a requirement 
on the beam duration T, with T > n/Q. Practically, this means the beam must rotate more 
than an half round, a condition which should be met obviously. With Q in the GHz range, 
this imply a beam pulse longer than a few nanoseconds. 

Besides this requirement, the quadrature f l47l ) is an oscillating integral vanishing for 
Q ^ oo. If it does so oscillating around 0, there will be an infinite number of values of 
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canceling the first harmonic. Consider the parabohc profile introduced by Eq. fllll) . The 
first harmonic amplitude is given setting n = 1 in Eq. ([I2D- Imposing = gives, 

_ = tan— . (49) 

These two functions are plotted on Figure M, and one can check the equation above has an 
infinite number of solutions. The first ones can be found numerically. Then, the solutions 
are found just before Xk = {2k + l)7r/2, with G N. Considering tanx ^ {xk — x)~^ near 
x/e, we can find an approximation for Xk- The "magic values" of the product QT/2 canceling 
the first harmonic are eventually given by, 

xi = 4.493, 
X2 = 7.725, 
xs = 10.904, 

By choosing any of these QT = 2xk^ the main asymmetry will come from the second har- 
monic only. The present procedure can be adapted to any beam profile and should allow to 
determine the optimum stability setup accordingly. 



VI. CONCLUSION 



We have conducted the Fourier analysis of the ion density deposited on the circular 
absorber by a rotating beam. For a point-like beam with temporal profile, calculation can 
be performed exactly, and the spectrum of the deposited ion density expressed in terms of the 
spectrum of the temporal profile. The extension of the calculation to a spatially extended 
beam can also be assessed exactly, and a geometrical form factor is found correcting the 
point-like spectrum. As long as the beam radius is smaller than the absorber radius, the 
ID spectrum, accounting only for the beam temporal profile, can be used for 2D beams. A 
stability analysis elaborating on the results from Piriz et al then gives the same result. 

It has also been proved that it is perfectly possible to cancel the first harmonic of the beam 
deposition around the absorber. This allows to further reduce the irradiation asymmetry, 
and render more realistic the use of a RTI planar interface analysis. Calculations have 
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been performed selecting a parabolic temporal beam profile, but they should eventually be 
conducted accounting for the beam profile chosen for the experiment. 
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Appendix A: Jacobian of Eq. (|22|) 

The Jacobian involved in Eq. (I22l) has to do with the change of variables {u^ 9^) to (i?, a). 
We start switching the 9 origin with respect to Fig. [5] so that 9 = (Rc,R). This does not 
bear any consequence on the present calculation because the differential elements remains 
the same by this simple transformation. Taking now the squares of the vectorial relations 
Rc + u = R and u = R — Rc gives, 



COS^,; = 



u = 



^R^ + Rl-2RRcC089 
R^ - Rl - 
2uRc 



(Al) 



Inserting the value of u from the first equation into the second one yields, 



cos ft 



R cos 9 — Rc 



(A2) 



U 



U 



The Jacobian matrix then reads 




(A3) 



The determinant of this matrix is found 



J{a, R) 



R 



(A4) 



^JR^ + Rl-2RR^cose 



where u{R, 6) as given by Eq. (lAll) is at the denominator, which proves 



u{R,e)J{R,e) = R. 
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(A5) 



This result may not seems surprising. We basically switch from one system of polar coordi- 
nates to another with different origin. However, the fact that in such case the Jacobian is 
so simple is not so well known, and we choose to include the proof. 
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